Negative Poisson's ratio materials via isotropic interactions 

Mikael C Rechtsman^, Frank H. Stillinger^, and Salvatore Torquato^'^'^'^ 
^Department of Physics, Princeton University, Princeton, New Jersey, 08544 

'^Department of Chemistry, Princeton University, Princeton, New Jersey, 08544 

00 ■ 

Q ' Program in Applied and Computational Mathematics 

O' 

CN ' and PRISM, Princeton, New Jersey, 08544 

"^ Princeton Center for Theoretical Physics, Princeton, New Jersey, 08544 <^i^d 

CN ! ^School of Natural Sciences, Institute for Advanced Study, Princeton, New Jersey, 08544 

(N 

(Dated: July 22, 2008) 

• 1— I'l 

c/3 ; Abstract 

We show that under tension, a classical many-body system with only isotropic pair interactions 
in a crystalline state can, counterintutively, have a negative Poisson's ratio, or auxetic behavior. We 
derive the conditions under which the triangular lattice in two dimensions and lattices with cubic 
symmetry in three dimensions exhibit a negative Poisson's ratio. In the former case, the simple 
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Q ■ Lennard-Jones potential can give rise to auxetic behavior. In the latter case, negative Poisson's 
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ratio can be exhibited even when the material is constrained to be elastically isotropic. 
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Materials with negative Poisson's ratio (NPR), the so-called "auxetics," are those that 
when stretched in a particular direction, expand in an orthogonal direction. NPR behavior 
is a counterintuitive material property that has been observed only in a handful of materials 
that often have intricate structures and characteristic lengths much larger than an atomic 
bond length. NPR materials have a great deal of technological potential, for example, to 
increase the sensitivity of piezoelectric transducers [l|; as components in microelectrome- 
chanical systems; and as shock absorbers and fasteners [2|, yj. 

A negative Poisson's ratio had not been observed in any elastically isotropic material until 
the discovery of certain foam structures with re-entrant structures Ml. NPR behavior has 
also been observed in multiple- length-scale laminate composites [5|, polymeric and metallic 
foams |6|, inverted honeycomb and other structures fabricated using soft lithography 2], 
and scaffold structures made out of springs, hinges, and rods ITJ]. It has been found in cubic 
atomic solids when they are stretched in the [110] direction |8|. 

In this Letter, we derive conditions for which NPR behavior is exhibited in classical 
many-body systems; this continues a research program in which inter-particle interactions 
are sought for targeted material properties. Examples of such inverse problems include 
optimization of pair potentials to give rise to negative thermal exp ansion [91], and a method 



to derive potentials that yield targeted classical ground states |10|. 

We report here that under tension, two- and three-dimensional systems with isotropic 
two-body interaction potentials can show NPR behavior in the crystal phase as long as cer- 
tain linear equalities and inequalities involving the interaction potential are satisfied. This is 
an unexpected result, since an inherently anisotropic behavior arises from isotropic interac- 
tions; indeed, most previously discovered NPR materials exhibit complex, carefully designed 
anisotropic interactions. We show this to be the case at zero temperature for the elastically 
isotropic triangular lattice in two dimensions, and for the fee lattice in three dimensions. 
In the latter case, NPR behavior is exhibited even when the material is constrained to be 
elastically isotropic. We first describe the calculation of the Poisson's ratio for any dimen- 
sion. Then, we present results for the two- and three-dimensional cases, including the elastic 
constants and NPR constraints. In order to demonstrate that NPR behavior is achievable at 
positive pressure, we present an example in which this is achieved by including three-body 
interactions in two or three dimensions. 

Consider a set of A^ -|- 1 particles, with positions {r„} [n > 0), that occupy a particular 



Bravais lattice in a state of zero strain. Under strain that is uniform throughout space, the 
new positions of the particles are 

x„ = (I + E)-r„, (1) 

where I is the unit tensor in d dimensions and E is the second-rank strain tensor with 
components Sij. The latter is constrained to be symmetric (i.e., Sij = Sji) in order to remove 
simple rotation. For simplicity, we take the origin to be at tq, and thus the energy per 
particle of the system can be written as 
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« = i^0(|x„|), (2) 

n=l 



where (^ is the pair interaction potential and A^ is the total number of particles excluding 
that at the origin. At zero temperature, the enthalpy and Gibbs free energy per particle are 
equivalent and equal to g = u + pv, where p = —du/dv, and v is the rf-dimensional specific 
volume. The equilibrium volume is found by minimizing g with respect to v at fixed p. In 
order to calculate the second-order elastic constants and Poisson's ratio, g is expanded to 
quadratic order in the strain tensor, Sij. The expansion is taken around equilibrium; thus, 
it contains no linear terms, and can be written as (7 ~ (^o + ^^ijki^ij^ki, where go is the 



zero-strain free energy, Xij^i are the second-order elastic constants of the system ll|], and 
the Einstein convention is used. The Poisson's ratio z/ = —et/sl is calculated by imposing 
a strain, el = Eijiiirij, and minimizing the free energy to find the strain in a transverse 
direction, et = Eijn[n'-, where n and n' are unit vectors in the original and transverse 
directions, respectively. 

The number of independent elastic constants is determined by the dimensionality and 
rotational symmetry of the lattice in question. For example, in two dimensions, square 
lattices have three independent elastic constants, and triangular lattices are "elastically 
isotropic" (i.e., elastic properties are independent of direction) and thus have only two |ll| . 
Lattices with cubic symmetry have three independent elastic constants, the least of any 
crystal structure in three dimensions. 

Two Dimensions: Triangular Lattice Example 



This analysis can be applied to any structure in two dimensions, but perhaps the simplest 
case is that of the triangular lattice. The free energy (neglecting qq, an irrelevant offset), is 
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where Ai and A2 are the elastic constants, and v is the two-dimensional "volume," namely, 
the area. Since both terms are quadratic invariants of the strain tensor, g is rotationally 
invariant. This property holds true of the Poisson's ratio, which is 
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The free energy expansion must be a positive-definite quadratic form, or else the lattice 
is unstable. Thus, the eigenvalues of the Hessian matrix of gn-i must be positive. This is 
true if and only if Ai > and A2 > 0. The Poisson's ratio must thus fall within the range 
— 1 < Vtri < 1- By expanding Eq. ([2]) as well as the area to quadratic order in the strain 
tensor, we obtain the elastic constants in terms of the first and second derivatives of the 
pair potential, evaluated at the neighbor distances of the lattice: 
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where the sum is over every particle except that at the origin, and Xi and yj's are the 
Cartesian coordinates of the particles in the unstrained lattice. A straightforward calculation 
shows that the Poisson's ratio, given in Eq. (jl]), can be written as 

1 + 2pK,T 
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where kt and p are defined by: 
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l7 = -ii:f=i|r#'(|r.|). 



Thus, the Poisson's ratio, given in Eq. i^, cannot be negative as long as p > 0, since for 
mechanical stability, kt > 0. Importantly, it is only for negative pressures, i.e., when the 
system is under tension, that utri can be negative. 

Even if the interaction potential extends only to the nearest neighbor in the triangular 
lattice and is zero beyond it, NPR behavior can be achieved at negative pressure, as long 
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FIG. 1: (Color online) Regime of negative Poisson's ratio in a triangular lattice, where the pair 
potential extends to the first neighbor only. The parameter space is composed of the first and second 
derivatives of the pair potential evaluated at the nearest-neighbor distance. In the continuously 
shaded region, vtri-, given in Eq. ([4]), is negative. In the grid-shaded region, the lattice is stable. 
The overlap region is given in Eq. ([8]). The pressure is positive to the left of the dotted line and 
negative to the right, 
as: 

0'(a) > 0, 0'(a) < act)" {a) < 50' (a) (8) 



where a = ^2^ /a/3 is the lattice constant. Figured] depicts these inequalities in a parameter 
space of 0'(a) and (j)"{a). In the case that the potential extends beyond the first neighbor, 
corresponding inequalities to those in Eq. ([8]) can be found, involving the first and second 
derivatives of the potential at relevant neighbor distances. 

To understand intuitively why NPR behavior only occurs at negative pressure, consider 
the function $(r^) = 0(r), which is the interaction potential rewritten in terms of the square 
of the inter-particle distance. A condition for lattice stability is $"(a^) > 0, which physically 
means that the "effective repulsive spring constant" between nearest-neighbor particles must 
decrease with distance. Thus, if the bonds between neighboring particles are thought of as 
springs with zero rest length, an imposed outward strain has the effect of weakening the 
effective spring constants of the bonds in the transverse direction. Under positive pressure, 
such a weakening causes a contraction, but under tension it causes an expansion. 

When is taken to be the well-known 12-6 Lennard- Jones (LJ) potential, given by 
4>Lj{f) = e [{b/ry^ — 2(6/r)^], there is a range of lattice constants for which NPR behavior 
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FIG. 2: (Color online) Region of lattice constants (indicated by the rectangular box) for which the 
Poisson's ratio is negative in a triangular lattice, using the LJ interaction potential, (pu- Pressure 
is positive to the left of the dotted line and negative to the right; thus, NPR behavior only occurs 
at negative pressure. To the right of the rectangular box, the lattice becomes unstable, 
is exhibited, namely 1.05966 < a < 1.08706, as demonstrated in Fig. [2l This corresponds to 
a pressure range of — 4.6363(e/6^) > p > — 4.8516(e/6^). Therein, the Poisson's ratio takes 
on values of through —1, the lower bound for stability. Note that a sufficient number of 
neighboring particles are included to accurately calculate the energy. That the Poisson's 
ratio is negative in this system is a surprising result, since the simple LJ potential has been 
studied extensively, and was not thought to exhibit this behavior. 

Other two-dimensional crystals have been studied for a general 0, including the square 
lattice and non-Bravais crystals such as the honeycomb and kagome. In all cases, NPR 
behavior was only found under tension. 

Three Dimensions: Lattices with Cubic Symmetry 



We proceed here with a calculation for lattices with cubic sjTiimetry that is analogous to 
that for the triangular lattice. The most common examples of these are the face-centered 
cubic (fee), the body-centered cubic (bcc) and the simple cubic (sc) lattices. Such systems 
have three independent elastic constants and are generally not elastically isotropic. The 
most general expression for the Gibbs free energy of cubic systems is. 
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where Ai, A2, and A3 are the elastic constants, and v is the specific volume. The Poisson's 
ratio is a function of both the direction of the imposed strain and the chosen transverse 
direction. As an example, if the strain is imposed along the [100] direction, the Poisson's 



ratio is 



''cubic — T r~r~- (10) 



In order for the Hessian of the quadratic form given in Eq. ([9]) to be positive definite, a 
necessary condition for lattice stabihty, we must have that Ai > 0, A3 > and — Ai/2 < 
A2 < Ai. If we expand Eq. ([2]) to quadratic order in the strain tensor components, we obtain 



Ai = -P+ ^Ef=i {^) [|r.|V'(|r.|) - |r#'(|r,|)] , 

A2 = f + ^ - i Ef=i {^Y [Ir^l V'(lr.l) - |r#'(|r.|)] , 

A3 = -P + ^Ef=i (a)' [|r.|V'(|r.|) - |r,|0'(|r,|)] , with 

^ = W. ^li (r-)' [|r.l V'(lr.l) - |r#'(|r.|)] , and 

P=-^Ef=i|r#'(|r,|). (11) 

where p is the pressure and kt is the compressibility. Note that these expressions apply 
equally well to any Bravais lattice with cubic symmetry. If we impose an additional linear 
constraint on the elastic constants, we find that we can impose elastic isotropy on the system, 
namely, if we enforce the following: 

Ai = A2 + 2A3, (12) 

then the system becomes elastically isotropic because the free energy given in Eq. ([9]) can 
be written as a function only of quadratic invariants of the strain tensor. Even with this 
additional applied constraint, we find that the system can exhibit NPR behavior under 
tension, with the Poisson's ratio in any direction given by Eq. flTOj) . The Poisson's ratio 
must fall between —1 and +1/2 in this case. If we consider a particular lattice we can, 
by employing Eqs. ( ITO|) . ( ITT1) and ( TT2|) find inequalities, which involve the pair potential 
evaluated at the neighbor distances, that describe the regime in which the Poisson's ratio is 
negative and the system is elastically isotropic. For each set of coordination shells included 
in the calculation, such inequalities can be found; none were found to allow for NPR behavior 
at positive pressure. 

Consider as an example the fee lattice, which has 12 nearest neighbors (at distance a/\/2) 
and 6 next-nearest neighbors at distance a, where a is the side length of the cubic cell. If the 
potential extends only to the nearest neighbor, NPR behavior is not possible. Both NPR 
behavior and elastic isotropy can be exhibited at negative pressure if the pair potential 



extends to the nearest and next-nearest neighbors, if the following constraints are satisfied: 



40'(a) - 4a0"(a) + a0" (^) = ^20' {^] 
^0' (^) < «0" (^) < 40'(a) + 5v^0' (^) , and 

40'(a) + 9v^0' (^) < 5a0" (^) . (13) 



At zero pressure, the Poisson's ratio goes to 1/4, as predicted by the Cauchy relations |12 |. 
Our analysis suggests that NPR behavior does not occur at positive pressures in crystals 
when the system contains only pair interactions, and the material is elastically isotropic. 
However, we present here a three-body potential that yields NPR behavior in close- packed 
two and three dimensional lattices by construction at zero temperature and positive pressure. 
In order to produce this behavior, the potential has a built-in energy cost associated with 
deforming the equilateral triangles in the two-dimensional triangular lattice and the three- 
dimensional close-packed lattices. The three-body potential is given by 

03(r, s, t) = af{r)f{s)f(t)F{r, s, t), (14) 

where a is a positive constant, r, s, and t are the side lengths of a triangle defined by a 
triplet of particles, / is some function that goes to zero sufficiently quickly that only nearest 
neighbors are within range of the potential, and 

(r + g + tf 

33/2 [2 (^2g2 + ^2^2 + ^2^2) _ ^4 _ ^4 _ ^4] J 

a function that is zero if the triplet of particles defines an equilateral triangle, but is greater 
than zero otherwise. Thus, if an outward strain is imposed on the system, then regardless of 
the ambient pressure, the system will expand in the transverse direction, if a is sufficiently 
large. 

In conclusion, we have shown that in two and three dimensions, classical systems with 
only pair potentials can have a negative Poisson's ratio at zero temperature, a surprising 



Fir, s, t) = — ^^^ ^ ^ - 1, (15) 



result. To the authors' knowledge, this has not previously been made explicit 13|. However, 



this auxetic behavior is only present when the system is at negative pressure, and thus not 
in thermal equilibrium. NPR materials may potentially be experimentally produced using 
synthetic techniques that rely on kinetic effects; examples include tempered glass |M], and 
even colloidal crystals [l3\. In two dimensions, it was proved that NPR behavior could 
be found at negative pressure; the proof was shown here for the triangular lattice, but a 



similar result also holds true for the square lattice. In three dimensions, a set of constraints 
on the pair interaction was found such that, if satisfied, the fee lattice is both elastically 
isotropic and has NPR behavior at negative pressure. Although the fee was chosen as an 
example, the calculation may be generalized given the expressions for the elastic constants 
of cubic systems reported here. We also presented a three-body interaction potential that by 
construction gives rise to a solid with elastic isotropy and NPR behavior at zero temperature 
and arbitrary pressure (negative and positive values). This suggests that the requirement 
of negative pressure is limited to systems with only pair interactions. A general proof of 
such a statement does not exist and will be considered in future work. In other future 
work, we hope to describe under what conditions NPR behavior can be observed in colloidal 
crystals under tension, given the experimentally realizable interaction potentials between 
colloidal particles. Finding NPR behavior over a wide range in temperature and pressure is 
a challenging optimization problem that we also intend to address. 
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